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1. Introduction 

Over the years a particular type of mathemat- 
ics has been developed for the evaluation of the 
integrals that occur in Feynman diagrams with 
loops. The solution of these integrals often in- 
volves the introduction of Feynman parameters 
and the manipulation of a special set of func- 
tions. Yet we are running into limitations im- 
posed by the complexity of the problem. Hence 
most quantities that arc important for comparing 
theoretical predictions and experimental observa- 
tions have not been evaluated beyond two loops. 
Yet it is important that a number of these quan- 
tities will be computed at the three loop level. 
Among these are the QCD structure functions in 
deep inelastic scattering. Considering the accu- 
rate measurements at HERA and the need for 
precise parton distribution functions at the LHC 
this may currently be considered as one of the 
most important calculations to be done. The con- 
ventional techniques have thus far not led to any 
results. Hence the introduction of different tech- 
niques may be called for. One method for the 
evaluation of structure functions has been around 
from the beginning of QCD [ Q[ ^, but has been 
used only very occasionally and then not to the 
fullest of its potential. This method is the eval- 
uation of the anomalous dimensions and/or the 
coefficient functions in McUin space in such a way 
that all even or all odd moments arc determined. 
When all even or all odd McUin moments of a 
function arc known the function of which these 
are the McUin moments is uniquely determined [ 
1). This is within a resonable set of assumptions 
about behaviour of the functions involved. Hence 
a systematic method to obtain these Mellin mo- 



NIKHEF-00-008 



ments would be sufficient to obtain the anoma- 
lous dimensions and/or the cocfhcicnt functions. 
This method has been applied several times in 
the past [ ^[ |^[ H, but the original computa- 
tion [ ^ of the full set of two loop coefhcient func- 
tions used a more conventional method. Here I 
will report on attempts to use this Mellin method 
even at the three loop level. This needs a very 
careful study of the functions involved. Much of 
this is new mathematics. We will address here 
harmonic sums, harmonic poly logarithms, Mellin 
transforms, inverse McUin transforms and differ- 
ence equations. 

2. Harmonic sums 

The basic functions that occur in Mellin space 
are the harmonic sums (HS) [ |l[ ||[ 0. These 
are defined by 
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while 'higher' sums can be defined recursively 
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The parameters m are called the indices, while 
N is called the argument. Wc define the weight 
of a sum by the sum of the absolute values of 
its indices. Whereas the two loop structure func- 
tions needed only sums of up to weight 4 [ pT|, 
we expect the three loop coefficient functions to 
need sums of weight 6. There exists an alterna- 
tive notation in which the indices take only the 
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values 1,0 and -1. In this notation a indicates 
that the nearest nonzero index to the right of this 
zero should be raised by one in its absolute value. 
Hence 
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In this notation one can see easily how many in- 
dependent sums there are for a given weight. All 
functions of a given weight w have exactly w in- 
dices. All these indices can have the values 1,0,-1 
with the exception of the last one, which can have 
only the values 1 and -1. Hence there are 2 3^""^ 
independent sums of weight w. One of the nice 
properties of these sums is that they form an al- 
gebra. This means that the product of two sums 
with the weights wi and W2 respectively can be 
expressed as a sum of terms, each with a single 
harmonic sum of weight w = wi + W2. 

These sums show themselves in a variety of 
ways. The most immediate is the expansion of 
the F-function. 
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At the positive side one has a slightly messier 
expansion. 

Another one is in sums of the type 
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In general we will need sums over these sums. 
This can be set up in a variety of ways. In some 
ways the multiple sums are so intertwined that 
we cannot resolve them. In other ways we can ar- 
range that we have always one parameter sums, 
meaning that for instance the upper limit of the 
summation and the summation parameter itself 
are the only parameters in the problem and hence 
the upper limit is the only parameter in the an- 
swer. If our sums can be expressed as nested 
one-parameter sums we make a good chance to 
be able to resolve them. Some examples of such 
single parameter sums that can be handled are [ 
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If we assume that each diagram can be expressed 
in terms of harmonic sums and we can bring a hi- 
erarchy amoung the diagrams in which the more 
complicated diagrams can be expressed as single 
sums over simpler diagrams, then the above con- 
dition is fulfilled. 

3. Harmonic polylogarithms 

In the end we would like to have functions of 
X = Q-Q /2P -Q rather than just the McUin mo- 
ments. When the two-loop structure functions 
are computed directly as functions of x, one runs 
into complicated integrals of which the solutions 
can be expressed in terms of dilogarithms and 
trilogarithms [ ^ of composite arguments. At 
the three loop level one could imagine that tetra- 
and penta-logarithms might do the job, but it 
turns out that these functions are not sufhcient. 



Also the Nielsen polylogarithms [ 13 1 are not suf- 
ficient. Hence we need the so called harmonic 
polylogarithms (HP) [ |lj] which are also closely 
related to multi dimensional polylogarithms [ |l^ . 
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in which rh^ stands for m^,, ■ • • , mi. One can de- 
rive a number of transformations for these func- 
tions that aUow a moderately easy evaluation of 
these functions in the entire complex plane with 
the exception of the points around +i and ~i 
where convergence is rather slow. We will need 
them on the interval — 1 only. When one makes 
a taylor expansion of these functions, the coeffi- 
cients are combinations of harmonic sums. 

One can define weights for these functions in 
the same way as was done for the harmonic sums. 
It shows that for each weight there are S"" differ- 
ent harmonic poly logarithms. Like the harmonic 
sums the harmonic polylogarithms form an alge- 
bra, although the rules are somewhat different. 
This makes it very easy to bring expressions into 
an unique standard form. 

4. Mellin transforms 

The Mellin transform is defined by 

M(/(x))= f dxx^'-'fix) (8) 
Jo 

It is not significant whether we use — 1 or 
iV — 2 in this definition. When we do such in- 
tegrals the result can be expressed in a combina- 
tion of harmonic sums in N and harmonic sums 
in infinity. The latter ones are just constants. 
Some efforts have been spent recently on reduc- 
ing these constants to a minimal set. This has 
been programmed up to weight 9 [ |l^ [ ^ , which 
should be sufficient even for the five loop anoma- 
lous dimensions. We will need them to weight 6 
in which case there are 'only' 8 different constants 
and combinations of them. 

The interesting part comes when we look 
at Mellin transforms of functions of the types 
H{x)/{l—x) and 7J(a;)/(l-|-a;). For a given weight 
m there are 3™ functions H and hence the above 
set of Mellin transforms has 2 3™ elements, which 
is identical to the number of harmonic sums of 
weight m + 1. Indeed there is a one to one cor- 
respondence between the above Mellin integrals 



and harmonic sums of weight m -\- 1. This is to 
say that each of these integrals results in a sin- 
gle harmonic sum of weight m + 1 and a number 
of harmonic sums of lower weight, possibly mul- 
tiplied by some of these constants. This allows 
for a method of inverse Mellin transformation, 
but there is still one little problem. When using 
the Mellin transform one obtains not only har- 
monic sums but there may also be a factor (— 1)^. 
This doubles the number of potential terms in 
the space of harmonic sums. Let us write it dif- 
ferently however. Instead of using the factors 1 
and (-1)^ we use the factors (1 ± (-l)^)/2. Now 
we only have to specify whether the expression we 
have is for even N of for odd N and there is a one 
to one correspondence again. This is fully in line 
with the theorem that knowledge of only the even 
moments or only the odd moments is sufficient to 
reconstruct the function of x. 

The method to do the inverse Mellin transform 
is now simple, provided one has a simple algo- 
rithm to say which index field of the harmonic 
polylogarithms corresponds to which index field 
of the harmonic sums. Such an algorithm exists [ 
0. 

1: Given a HS, construct the single term with a 
HP with the highest weight and the denominator 
1/(1 ± x). Subtract and add. 

2: Do a Mellin transform on the subtracted 
term. This will create a term that cancels the 
original HS and gives also terms with lower weight 
HS's. 

3: Keep repeating this procedure untill there 
are no more terms with HS's. The remaining con- 
stant terms will get a factor 5{1 — x). 

5. Difference equations 

When we study diagrams as a function of iV, 
usually one cannot come up with a simple expres- 
sion in a direct constructive way. Instead one can 
construct equations of the type 

An{N)F{N) + Ai{N)F{N-l) + 

■■■ + A^{N)F{N^n) + G{N) = (9) 

Such an equation we call an n-th order difference 
equation. A simple example would be a first order 
equation which is also called a recursion. In that 
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case the solution can be written down as a single 
sum: 
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Under the right conditions the products can be 
written as combinations of F- functions in N,i,j 
and e. 

Recently other authors have also obtained dif- 
ference equations for the evaluation of diagrams. 
Tarasov [ |l^ obtains equations in the dimension 
of space-time. Laporta [ |l^ has a very promis- 
ing and systematic approach that has similarities 
with the above method, except for that it does not 
involve the Mellin parameter N. This results in 
a completely different approach to solving them. 

When we define basic building blocks (BBB) 
as diagrams in which the momentum P of the 
parton flows only through a single line in the di- 
agram, we have to worry about two BBB's when 
we do a two loop calculation and 17 BBB's when 
we do a three loop calculation. Of these 17 seven 
are two loop diagrams with a one loop insertion. 
The two BBB's at the two loop level result in 
first order difference equations. This was already 
noticed in [ ^. The 17 BBB's at the three loop 
level split into 4 first order equations, 12 second 
order equations and one third order equation. We 
have been able to set up all these equations. It 
is possible to define a hierarchy amoung the in- 
tegrals in which a more complicated diagram F 
is related in this way to simpler diagrams inside 
the function G. For example the equation for a 
ladder diagram may have two loop diagrams with 
a one loop insertion inside the function G. The 
coefficients Ai are usually polynomials in N and 
e = 2 — D/2. Let us assume that the function G 
can be expressed in terms of harmonic sums. For 
a number of types of these equations one can then 
make a constructive solution for the function F. 
This is usually rather complicated and moreover, 
we did not find constructive methods for all the 
equations that we derived. There does however 
exist a quite good method that does the job for 
all equations. This method goes as follows: 



• Divide the equation by N"^ in which q is the 
highest power of N occurring in the coeffi- 
cients Ai. 

• Synchronize the arguments of the sums in 
G and the denominators. 

• Write S„-,{N)/NP = Sp^^iN)- Sp^N -I) 
and similar forms when (—1)''^ is involved. 

• Assume that F(A^) —J2^i,n,m<^^S^{N-n). 
We have to include powers of ^3, C4 and (5 
as well, but let us ignore that here. We have 
to guess a bit what range of values we need 
for i and n. 

• We also bring the Aj {N)F{N—j) terms into 
the standard form with just a single S'-sum 
without denominators. 

• This defines a large number of equations in 
the Ci^n,rn which we solve. If the range in 
i and n was large enough we should find 
a solution. The maximum weight we use 
defines (amoung others) to what power in e 
our answer can be correct. 

• There will still be a few free parameters left, 
as we have not fixed the boundary condi- 
tions yet. For an n-th order equation we 
will need n values. These can be computed 
with the Mincer program, as they concern 
fixed values of N. 

• The boundary conditions give a few more 
equations, and these fix all remaining con- 
stants. This then determines the function 
F. 

This method is rather general and relatively fast. 
For the cases for which we could compare it 
with the constructive method it was significantly 
faster. It allows us to solve all equations we run 
into provided we do not run out of computer re- 
sources (in our method of solution there is a lim- 
itation in the 32 bits version of FORM [ |9| that 
limits the equations to between 1000 and 2000 
terms each. The 64 bit version has no realistic 
restriction in this respect. It would also be possi- 
ble to choose a less efficient way of solving these 
equations if the need were to arise.) 



5 



Because we have to solve these equations aU the 
way to weight 6 sums (of which there are 486 for 
each n) and for several values of n, these equa- 
tions involve often a few thousand coefRcients 
to be solved. Fortunately the system is rather 
sparse. 



6. Status 

We have used the equation solving program al- 
ready to solve all equations that occur at the two 
loop level in which the two loop integral has all 
lines with integer powers, but one line has a power 
involving e as is to be expected from three loop 
integrals of which one loop has been integrated al- 
ready. In addition we have constructed all equa- 
tions for the basic building blocks at the three 
loop level. Because it is still much work to pro- 
gram the whole hierarchy we are still far away 
from completing the computation of the simplest 
function G, solving the next function F, com- 
puting the next G, then the next F, etc. all 
the way to the most complicated diagrams. The 
most complicated nonplanar basic building block 
is determined by a third order difference equa- 
tion, while the most complicated benz type dia- 
gram has a second order equation that contains 
about 200000 terms before the integrals in the G- 
function are evaluated. After evaluation it should 
simplify considerably of course. There are how- 
ever no more fundamental problems left to obtain 
solutions for the basic building blocks. 

The reduction for composite diagrams into ba- 
sic building blocks follows a similar scheme. Take 
for instance the diagram 
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We derived a general formula for such diagrams 



that is relatively simple: 
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This allows us to express the NOi^ diagram into 
a simpler diagram by means of a simple recur- 
sion. Hence this can be solved, once we know 
the simpler diagram. We could either use the 
general program for difference equations for this, 
or write the solution as a single sum over this 
simpler diagram. This sum involves F-functions 
which have to be expanded (and give harmonic 
sums again) and in the end we can use some gen- 
eral procedures for such sums. The main problem 
is currently to derive all the equations and recur- 
sions that break all composite diagrams down into 
basic building blocks. This is just an enormous 
amount of work. After that all relations have to 
be programmed and tested which is another enor- 
mous amount of work. 
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